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Hughes & Gaztanaga (2001, see article in these proceedings) have presented reahs- 
tic simulations to address key issues confronting existing and forthcoming submm 
surveys. An important aspect illustrated by the simulations is the effect induced 
on the counts by the sampling variance of the large-scale galaxy clustering. We 
find factors of up to ~ 2 — 4 variation (from the mean) in the extracted counts from 
deep surveys identicpJ in area (~ 6arcmin-^) to the SCUBA surveys of the Hubble 
Deep Fields (HDF) Here we present a recipe to model the expected degree of 
clustering as a function of sample area and redshift. 



1 A model for the angular clustering 

Fig.|l| shows fluctuations in the galaxy counts as measured in the APM Galaxy 
Survey lI. Symbols with errorbars show the square root of the variance W2{d), 
e.g. AN/N = \/ui2, measured in squared cells of area A = 9^. The data can 
be described by a power law: 

AN/N = ~ {A/Aof 

with Ao ~ 7.6 X 10"^ sq.deg. and f3 ~ -0.175, for A < 2.7 sq. deg. For larger 
areas there is an exponential cut-off with a characteristic scale of Ac ^ 110 
sq. deg. The complete model is represented by: 

^ ^ (A/Aof exp[-AMJ, (1) 
and is shown as a solid curve in Fig.fl. 



2 Shot-noise 

The effects of shot-noise on the counts, 6.17. due to the small number of 
sources in the sample, is easy to quantify. If the mean number density of 
sources at a given flux is Af, then the number of sources in a sample of area 
A is: N — AJ\f. The total variance on such area is: 

U,tfal^-^^^ (2) 
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Figure 1. Mean square root deviation in number counts as a function of the sampled area 
in the sky. Circles with error-bars correspond to counts in square cells of different areas 
measured in the APM survey. The solid curve shows the model in Eq.[Q. The predicted 
fluctuations, scaled to a depth of X> ~ 2500/i~^Mpc {z ~ 3), are shown for different cluster- 
ing evolution models: fixed in co-moving coordinates (dotted line), stable clustering (dashed 
lines) and linear theory (long-dashed lines). The dot-dashed line includes the shot-noise 
contribution for a surface density of 2000 sources/sq. deg. 



where W2 is the intrinsic variance (in the case of high density) 0. Thus, the 
number counts variations due to intrinsic clustering and shot-noise is: 



AN 



1/2 



(3) 



where W2[ 9;z;T>\ is the angular variance of depth T) at redshift z. We next 
need to quantify the effects of projection and clustering evolution. 
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3 Projection effects 



Let us first assume that the 2-pt function ^2(?'i2) does not evolve in co- moving 
coordinates: £,2{ri2) — £,2{xi2), where ri2 = 2:12/(1 + z). We then have that: 

W2{0) = / dxidx2 ipixi) i;{x2) 6(3=12), (4) 

where dxi is the co-moving volume, ip{x) is the normalized probability that 
a galaxy at a coordinate x is included in the catalogue and 1^ is a cone (or 
pyramid) of radius (or side) 9 and infinite depth (with solid angle equal to the 
sample area A). In the small angle approximation, e.g. when the transverse 
distances are much smaller than the radial ones, we can relate the angular clus- 
tering of samples projected at two different characteristic co-moving depths 
T>2 and Vi by: 

W2[e ; Vi] = ^ iU2[ eV2/Vi ■V2]. (5) 

We take V2 ~ 400/i~^Mpc as the mean APM depth and Vi ~ 2500/i-^Mpc 
as the mean co-moving depth of a typical submm sample (the actual num- 
ber depends on both the mean redshift and the cosmological parameters, as 
given by the luminosity distance relation). Thus the curve that fits the APM 
measurements (solid curve) in Fig.|l| must be moved left by a factor 7)2 /T>i 
to account for the fact that, at larger radial distances, the same physical 
length subtends a smaller angle. The APM curve must also be moved down 
by a factor 7)2 /T^i because more galaxies are seen in projection at greater 
radial distances. The resulting prediction in this case, scaled from Eq.j|] to 
Vi ~ 2500/i^iMpc with Eq.[|), is shown as the dotted line in Fig.|l|. The 
dotted-dashed line includes the shot-noise contribution in Eq. for a submm 
850/Ltm flux limit of 3 mJy, with surface density A/" = 2000/sq. deg., which are 
characteristic of the HDF counts. Note how in this case the variance of the 
counts are dominated by shot-noise at almost all scales. These predictions 
agree well with the source-count variations found in our submm simulations 
of the HDF (see Hughes & Gaztafiaga in these proceedings), which also have 
clustering fixed in co-moving coordinates and the same mean depth. 

4 The redshift evolution of clustering 

We next model the redshift evolution of ^ in proper coordinates ri2 = 
a;i2/(l + z) as: 

e(ri2;z) ^ (l-Hz)-(3+^' C(^i2;0) (6) 
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For stable clustering (pattern fixed at proper separations) we have e ~ at 
small scales 13. For a power-law correlation with slope 7 ~ 1.7, linear theory 
gives £ = 7— 1~0.7 and non- linear growth gives e ~ 1. If the clustering 
is fixed in co- moving coordinates, then e ~ 7 — 3 ~ —1.3 which produces 
even less evolution than in the stable clustering regime. This fixed clustering 
model describes galaxies which are identified with high density-peaks: peaks 
move less than particles, which results in less evolution. Weak evolution is 
consistent with the strong clustering observed in Lyman-break galaxies at 
2 ~ 3, which is comparable to the clustering of present-day galaxies B. The 
above model of cluster evolution projects as: 

W2[e;z,;V]^ W2[9;z2;V] (7) 

where 7 varies from ~ 1.6 over the smaller scales to 7 ~ 2 near the exponential 
break. Here we use Z2 — 0.15 for the APM depth and zi ~ 2 for the depth 
in sub-mm samples. The resulting predictions, scaled from Eq.[0 to Vi ~ 
2500/1""'^ Mpc, are shown as the dotted (e = 7 — 3), short-dashed line (e = 0) 
and long-dashed lines (e = 7 — 1) in Fig.|l|. 

5 Discussion 

In summary, our recipe for the angular variance of number count fluctuations 
(t^)^ = W2 in a galaxy sample of area A = 6^, mean co- moving depth V and 
mean redshift z, is given by: 



where /3~ -0.175, Oc ~ 10.5 deg., 6*0 ~ 8.7 x 10"^ deg., 3 + e-7is between 
and 2, depending on clustering evolution, and J\f is the mean number density of 
sources at the given flux that produces shot-noise fluctuations (e.g. Eq.j^]). 
Several cases for the above model are shown in Figj^. Assuming the shot- 
noise is negligible, e.g. Af 00, to reach a 1% level of fluctuation in N 
(lower horizontal line in the figure) we require a sample of about A ~ 6 sq. 
.deg. for the model with co-moving evolution and about A ~ 1 sq. deg. 
for a inodel with strong clustering evolution. In the submm survey of the 
HDF □, with A = 0.001 sq. deg. (the left vertical line in the Fig. 1), we 
expect AN/N ~ 0.2 in the case of weak clustering evolution and ~ 0.05 for 
the strong evolution case. Thus, in principle, comparing number counts in a 
few more submm surveys, similar in area to the HDF, would provide a clear 
discrimination between clustering evolution models. 
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Nevertheless in small submm surveys the number counts for bright sources 
will be low, e.g. N = (PM ~ 3 for S'sso/jm — 3mJy in the HDF survey 
and hence the shot-noise correction dominates and masks any evolution in the 
clustering. This situation is shown as a dot-dashed line in Fig.|l| for the submm 
HDF survey. In this case, we could still subtract the shot-noise contribution 
using Eq. ^ although this will introduce large uncertainties. Thus, in the 
future, it is necessary to extend the surveys to larger areas and/or to lower 
flux densities in order to have a discriminating measurement nf clustering. 
The Gran Telescopic Milimetrico/Large Millimeter Telescope □ will be the 
optimal facility to provide this new generation of deep and wide mm surveys. 
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